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Abstract 

We have developed a theory of polymer entanglement using an extended 
Cahn-Hilliard functional, with two extra terms. One is a nonlocal attrac- 
tive term, operating over mesoscales, which is interpreted as giving rise to 
entanglement, and the other a local repulsive term indicative of excluded 
volume interactions. We show how such a functional can be derived using 
notions from gauge theory. We go beyond the Gaussian approximation, to 
the one-loop level, to show that the system exhibits a crossover to a state 
of entanglement as the average chain length between points of entanglement 
decreases. This crossover is marked by critical slowing down, as the effective 
diffusion constant goes to zero. We have also computed the tensile modu- 
lus of the system, and we find a corresponding crossover to a regime of high 
modulus. 
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I. INTRODUCTION 



While it has long been known that entanglement in homopolymers has an important 
effect on its strength, a thoroughly satisfactory theory of polymer entanglement is still a 
topic of current research. The classic experimental wort of Moore and WatsonS showed that 
the bulk modulus of cross-linked natural rubbers depends inversely on the average chain 
length (iV c ) between cross-links in the system, and that end corrections become negligible as 
the total average molecular weight /i gets very large. They pointed out the analogy between 
chemical cross-linking and physical entanglement. Thus their work applies in a qualitative 
sense to entangled systems as well. Their work extended the earlier pioneering work of Flory 
et all. They showed that at 100% extension, the modulus for vulcanized rubber was higher 
than the one calculated from Kuhn'sil result, obtained from Gaussian chain statistics. 

Edwards developed the tube theory of the effect of entanglement on elastic modulii of 
homopolymers using deGenne's idea of reptationi. This theory showed how entanglement 
enhances the tensile modulus of a homopolymer. He also developed a more detailed model 
of entangled ring polymers using notions from knot theory!. The basic idea behind this 
theory is an analogy between certain mathematical invariants describing intertwined loops 
and magnetic fields induced in wires by current-carrying loops. Prager and Frischi worked 
on this notion as well, as did Koniaris and MuthukumaiH 

More recently, interest has turned towards computer simulations of polymer networks, 
involving various levels of molecular detail, to understand the effect of entanglement on the 
strength of homopolymers. As examples, we mention the work of Termoniai and Biceranoi, 
who use phenomenological models of polymer networks to study their viscoelastic properties. 
Comparison with experimental data shows a varying degree of success, depending on the 
particular system studied. Holtzl et all use the more basic fluctuating bond theory to model a 
network of polyethylene strands to show that entanglement leads to non-afnne displacements 
under large tensile strains. 

In an earlier paper,0 we developed a gauge theory of self-assembly, and utilized renor- 
malization group ideas to study the onset of self-assembly in diblock copolymers. In this 
paper, we shall pursue a similar continuum approach to understand entanglement. 

Intuitively, one can see that entanglement could be described by assuming two extra 
terms in the Cahn-Hilliard functional,!^ one of which is a nonlocal attractive term, which 
gives rise to entanglement and the other, a soft-core local repulsive term which arises from 
the fact that the strands cannot cut across each other. We connect the parameters which 
appear in our theory to the underlying chain parameters with a simple model. We have 
shown (see Appendix A) how such a functional can be derived naturally using notions from 
gauge theory. 

The results derived from our continuum formulation will be seen to be reminiscent of the 
chain-theory approaches of Kavassalis and NoolandiB for flexible polymer networks, and 
that of Kroy and Frejil! for semi-flexible networks. They utilized a mean-field approach 
to locate the transition to the state of entanglement. Our theory is also somewhat similar 
to the paper of Castillo and Geldartlll, who use a (jr field theory, coupled to the replica 
trick (in the mean field approximation) of Deam and Edwards^! to study the vulcanization 
transition. 

We shall utilize a field theoretic approach in this paper to show that the onset of the state 
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of entanglement is a crossover phenomenon, rather than a pure phase transition, in that the 
effective diffusion constant goes to zero at the transition point, but the correlation length 
and the structure factor do not diverge. We shall go beyond the Gaussian approximation to 
locate the critical chain length at which the transition to the state of entanglement sets in. 
A physical reason that fluctuations become important near the onset of the state of entan- 
glement is that the average chain length between points of entanglement gets smaller. This 
underscores a difference between vulcanization and entanglement- An entangled network of 
polymers is more dynamic than a vulcanized network. The mean field approximation is 
expected to be correctS for the vulcanization transition. 

We have also computed the tensile modulus of the system. Corresponding to the critical 
slowing down discussed above, we also find a crossover in the modulus to a regime of high 
values. Limitations of the approximations made are discussed. 



II. A FIELD THEORY OF ENTANGLEMENT 

The continuum, mesoscale approach adopted in this paper assumes that we have per- 
formed some spatial averaging of our polymeric system, so that the order parameter is the 
local concentration of the polymers. Our mesoscopic theory of entanglement in polymers 
is based on the intuitive notion that physical entanglement can be captured by a non-local 
attraction between the polymers which causes them to remain in proximity. There must be a 
balancing, repulsive local energy term which says that the polymers cannot cut across each 
other. The starting point of our mesoscale theory is an internal energy functional which 
is quadratic in the gradient of the local number concentration. For the moment, we will 
consider isolated systems, so that the quantity that is conserved is the internal energy^. 
We will shortly consider entropy effects as well. Consider the following form for the energy 
functional: 

PU = 13 ! u (c(s))d 3 s (1) 



**<■»- (§)^ < 3 » 

where repeated indices are summed over, s is a dimensionless co-ordinate variable, k is 
Boltzmann's constant, T is the temperature, and c is the number concentration of the specie. 
The local concentration c is normalized by dividing by some characteristic inverse volume. 
The constant g is analogous to a dimensionless diffusion constant. Such energy functionals 
have been considered over many years as contributing to the total internal energy of both 
unary and binary mixtures.0 We will use this form as our starting point to suggest a more 
complete energy functional: 

(3U eS = PU + (£) I d 3 s c( S )c(s) - (2. )U I dV c(s) eXpM|s ~ s/|) c(s') (4) 



2 I J ww y 2n j j j w , _ g , 
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where a 2 , 7, 5 are positive constants. The local repulsive term is indicative of the fact 
that polymers cannot cut across each other. This is in effect a soft-core repulsion term, and 
the softness arises because we are studying a homopolymer network at a mesoscale, where 
polymers may pass by each other, without actually cutting across each other. The nonlocal 
attractive term gives rise to entanglement, as it causes portions of the network within the 
screening distance 1/5 to be attracted to each other. Equation ^ is the basic statement of our 
theory. Note that the two terms we just discussed have signs opposite those of corresponding 
terms in theories of self-assembly0. 

In what follows, we shall set 7 = a 4 , and 5 2 = \^2a, with a 2 = g 2 /2. A strong 
motivation for this choice of parameters is provided in Appendix A, where we use notions 
from gauge theory to derive Eqn.f|, with the parameters having the forms given above. 
Another explanation for such a choice is as follows. With our choices for the parameters, 
U e ff in momentum space may be written as: 



U eS = J -0-c*(k)[-V2ak 2 - V2ak 2 /(1 + V2k 2 /a)]c(k) (5) 



where the carats indicate a Fourier transform. Thus we see that the choices made for 
the parameters are equivalent to generalizing the diffusion constant g = \^2a — > v2at[l + 
1/(1 + \^2k 2 /a)], i.e., a non-local diffusion constant is obtained. If we now extremize the 
functional, the Euler-Lagrange equations can be written in conservative form as: 

V- /(S)=0 

m= V / -^exp(ik-s)[l + 1/(1 + V2k 2 /a)]c(k) (6) 

where I(s) can be interpreted in the conventional manner as a mass current. The 
divergence-free nature of this current makes it clear that with our choice of parame- 
ters, our internal energy functional preserves number conservation. This is quite ap- 
propriate, since the internal energy is the quantity which is conserved for isolated sys- 
tems. For an arbitrary choice of parameters, the Euler-Lagrange equations have the form: 
V • I'(s) = Source/Sink Terms, indicating that number conservation can be a problem. 

While our choice of parameters may appear to be overly restrictive, it turns out to be 
sufficiently rich to provide a description of the onset of entanglement in polymers. We will 
not explore more general sets of parameters in this paper. 

Before we can compare our theory with experimental data, we need to consider the fact 
that our system is not really isolated, and may be in contact with an energy reservoir, 
perhaps as it is being acted on by mechanical forces in a stress experiment. For a system 
in contact with an energy reservoir, the quantity that is conserved is the Helmholtz free 
energy^ A = U c s — ST, where 5* is the entropy of the system. In general, when a system is 
in contact with an external reservoir number conservation is not always be gauranteed. As 
an example, we point out that in stress-strain experiments, a polymeric sample is clamped at 
two ends in such a way that individual polymers may leave the sample volume being tested. 
For this reason, we will not impose number conservation on the Helmholtz Free energy, as 
we did on the internal energy U e g. The entropy of our system will be written in the usual 
form: 
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d 3 s c(s) ln[c(s)] (7) 



This entropy term provides the free energy a single minimum. To ease computations, we 
shall expand cln(c) in a power series about the characteristic inverse volume Z~ 3 (= 1 in our 
dimensionless units), retaining terms up to fourth order: 

2 3 4 

(l + c)ln(l + c)^c+^-^- + ^ (8) 

This expansion yields one minimum, just as the exact expression above for the entropy 
(Eqn. |7|). Consequently, we do not expect this system to display a phase transition, but 
rather a crossover from an un-entangled state to a state of entanglement. Finally, we note 
that in our present theory, entropy yields the crucial nonlinear terms which will describe 
the crossover to a state of entanglement, in contrast to our gauge theory of self-assembly0 
where entropy did not play a dominant role. 

We define the two-point Green's function as usual via S(x, x') = (5 2 / ' 8J(x)5J(x'))Q[J}, 
where Q[J] = JVc9(l + c) exp[-/3(C/ eff - ST) - f d 3 sJ(s)c(s)}, where 9(1 + c) is a step 
function that indicates a restriction to physically acceptable values of the concentration. 
In practice, we shall be restricting our attention to small deviations of c from its averge, 
so that the step function is implicitly accounted for during calculations. In the quadratic 
approximation, the structure factor is: 

^ = 1 + a'k 2 + a' 2 fc 2 /(l + 2k 2 /a') (9) 

where a' = \/2a. Equation [5] displays a peak at the origin, as one might expect from the 
fact that entanglement creates blobs which are distributed at random within the system. 
The width of the peak indicates an inverse of the correlation length between the blobs. 
With this physical interpretation, y/a is a measure of the distance between concentration 
fluctuations (i.e. between points of entanglement). The decay of S(k) is affected by the 
value of a. As a decreases, the structure factor looks more diffuse. Thus, a decrease in 
a signifies a shift to a state of higher entanglement, as the concentration of entanglement 
points increases. 

Our results can be understood compactly in terms of the parameter a. This parameter 
is related to the properties of the underlying chains in the system. We offer below a simple 
model which connects a with some of the chain parameters. We offer this model as an 
example, in order to understand the basic physics in our theory. More sophisticated choices 
may be required for specific systems. We pointed out that g > plays the role of a diffusion 
constant. We expect g to decrease as entanglement increases in the sytem. Let N be the 
average number of links in the system. It is tempting to consider g = \/2kN~ 2 following 
deGennesB scaling arguments for the diffusion constant, k is a dimensionless constant. 
The parameter a ~ N~~ 2 , so that a decrease in a signifies a shift to a regime of higher 
entanglement, since the larger the average molecular weight of a polymer, the smaller is 
the value of N e . It must be emphasized that a ~ iV -2 is but a simple model, and that 
in general a = a(N) such that a(N) decreases as iV increases. It may be possible to 
obtain this dependence of a by comparison with experiments on specific polymers. Let us 



5 



now make some more definitions, viz., N e is the average chain length between consecutive 
points of entanglement, the entangled chain number density c e = pN avogadro / (p, N e ) , and 
the monomer number density c = pN avogadro / p , where p is the mass density of polymer, 
Navogadro is Avogadro's number, and p is the molecular weight of the monomer. I is the 
length scale in our theory and we shall take it to be Z = Ac" 1 / 3 , where A is a parameter which 
could be used to improve agreement with experiment. Such a phenomenological approach 
is analogous to the one in applications of reptation theory, where the tube diameter is 
often used as an adjustable parameter!. In what follows, we shall simply use A = 2 1 / 3 for 
exposition, as it leads to an expression for the tensile modulus in the Gaussian approximation 
which agrees with Kuhn's result. We shall use the length / to scale all other lengths in the 
system. These ideas are slightly similar to Stillinger's in another context0S In this manner 
we have attempted to relate our theory in an intimate fashion to the notion of entanglement. 



III. BEYOND THE GAUSSIAN APPROXIMATION. 

We shall now use diagrammatic methods to go beyond the Gaussian approximation to 
the structure function described at the end of the previous section. The reason is to be able 
to describe the crossover to a state of entanglement. As discussed in the previous section, 
the onset of entanglement is not a phase transition, but simply a crossover. 

Figures la and lb show the basic vertices in our theory. The figure captions describe the 
Feynman rules which go with these vertices. We shall compute only the first non-vanishing 
terms which arise from each of these vertices. The first order contribution of the cubic 
term is zero, as follows from symmetry considerations. We have to go to the second order 
in the cubic term to obtain a tadpole diagram which is non-vanishing, as shown in Figure 
2a. It serves to renormalize the correlation function in the long wavelength limit. Figure 
2b is the other setting sun diagram which comes from the second order contribution of the 
cubic term. It may be expanded in powers of its argument k. The term proportional to k 2 
helps to renormalize the diffusion constant g, and seves to diminish it, as one would expect 
entanglement to. Figure 3 shows the conventional bubble diagram coming from first order 
perturbation theory with the quartic term. It serves to renormalize the correlation function 
in the long wavelength limit. 

In order to render the integrals in our theory finite in three dimensions, we shall use 
the following regularization scheme. We shall perform an expansion of the denominator of 
the Gaussian structure factor in powers of k. We shall retain terms upto 0(k 6 ). This is 
essentially an expansion in inverse powers of a. This expansion yields the requisite higher 
order terms in the denominators of the Green's function to render our integrals finite, while 
ensuring that S(k) > 0. This method has the advantage of retaining the correct long- 
wavelength behavior, at the expense of high momentum behavior. This is acceptable, since 
we do not expect our theory to be correct at small wavelengths in any event. Our single- 
particle Green's function in the Gaussian approximation is now taken to be: 

^ = l + 2a'fc 2 -2A; 4 + /c 6 /(2a') ^ 
With this definition, the net contribution from diagrams shown in Figures 2a and 3 is: 
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S 2 « +3 («) = - ( 7 ) ^(O) (11) 



,4 

Figure 2b yields a /c-dependent contribution to the self energy: 



l\ r d 3 k' 

Sa + 5gk 2 + 0(k 4 ) (12) 



where: 



Sa(a) 



128 2 1 /%a 3 / 2 
/ 1 \ / 5 



(13) 



where the integrals were performed by approximating the denominator of the Gaussian 
Green's function by terms upto 0(k 2 ), as this suffices to guarantee convergence of the in- 
tegrals, so that there is no sensitivity to the higher order terms neglected. With these 
expressions, we see that the renormalized value gn = g — 5g of the diffusion constant de- 
creases as a is decreased, qr is zero near a = 0.18. Note that a decreases as we decrease 
N e the average chain length between consecutive points of entanglement. We thus see that 
as entanglement increases, the effective diffusion constant decreases, analogous to critical 
slowing down. This effect is similar to the considerations of Broderix et alj^l, who study 
the vulcanization transition in the mean field approximation, and find the diffusion constant 
going to zero as the vulcanization transition is approached. They point out that this re- 
sult agrees with experimental results. Given the analogy between vulcanization (chemical 
cross-linking) and physical entanglement, we believe this result should be experimentally 
observable in entangled sytems as well. Broderix et al obtained g^ — > linearly with the 
average concentration. We have obtained a more complicated dependence on the concentra- 
tion. We find that a ~ 0.18 when the renormalized diffusion constant goes to zero, yielding 



a critical N* « Jk/0.18. Since there is experimental evidence that entanglement sets in 
at N ~ (9(10 4 ), we deduce that k ~ O(10 7 ). Alternatively, one could begin by estimating 
K using results from the next section on the tensile modulus of polymers and experimental 
values for polymeric elastic modulii, and then obtaining a value for iV*. 

The origin of 5g > can be traced back to the nonlocal attractive term in U c g, defined in 
Eqn.[|. This nonlocal attractive term, which we interpreted as giving rise to entanglement, 
is responsible for a physical signature of the onset of entanglement, with g^ — ► 0. 



IV. TENSILE MODULUS. 



It is well-known that Kuhn's result for the tensile modulus, while yielding the correct 
trend, does not agree with experimental dataE on modulii by a large factor. Edwards' 
application of deGenne's reptation modeP provides an enhancement factor over Kuhn's 
result, and shows conceptually how entanglement leads to an increase in the stiffness of 
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the homopolymer system. We will show in this section how to obtain a similar result in 
our continuum treatment. More importantly, we will show we can go further, and describe 
a crossover, as the mean chain length between entanglements is decreased, to a regime 
where the tensile modulus, instead of remaining fairly constant, begins to increase extremely 
rapidly as a function of decreasing N e . The reptation model is unable to accomplish thisp 
as it assumes that the system is already in the entangled state, and does not account for 
inter-chain interactions, beyond assuming a preformed tube. 

The Helmholtz free energy in the Gaussian approximation is given by MS 

jr G = -\kTVc e S{k = 0) = -\kTVc e f d 3 xS(x) (14) 

We may represent a strained state of the sytem by the transformation x — > x! = x + u(x) 
in the above equation. This is possible because S(x) in the above equation represents the 
density-density correlation function < c(f)c(f— x) >, so that when the system is strained, 
c(r — x) in our theory shifts to c{f — x'), where f = f + u(r). For the case of homogeneous 
deformation, we shall take u(x) —*e-x, where e is the strain tensor. The strain is assumed 
to be volume-preserving, so that d 3 x = d 3 x'. Our approach is similar to that of Castillo and 
GoldbartEl. It is now easy to show, using a Taylor series expansion that the change in the 
pressure P = — (dJ r /dV) T N is: 

AP G ^ ~kTc e S(k = 0)e aP e rS V a ^s + 0(T) 

/ D af 3~ / s= (S a pSjS + S ay Sf3s) (15) 

where the subscript G denotes the Gaussian approximation. As needed, we can consider 
the expansion of Free energy to include higher orders of the strain tensor.^1 APq is a measure 
of the change per unit volume of the energy of the system under strain. In analogy with a 
simple harmonic oscillator, the force tensor which constrains the system from undergoing a 

strain *e is given by —kTc e S(k = 0)*e : T>. Thus, the stress a required to produce this strain 
is: 

<-> 

a = kTc e S{k = 0)~ : V (16) 

One may now readily write down the tensile modulus as: 

Yg = kTc e (17) 

This is identical to the well-known result obtained by Kuhn. And it's origin is purely 
entropic. Our goal is to go beyond Kuhn's result, and to do that, we shall dress the bare 
propagator <S(k) using the diagrams shown in Figures 2 and 3. This immediately leads to 
the renormalized result Yr\ 



Y R (a)= kTc e S R (k = 0) 

1 



kTcp 



1 - S 2a +3(") - 8a(oi) 



(18) 
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The first of these equations is similar to the connection made between the structure 
factor in the long wavelength limit and the bulk modulus by Kirkwood@. 

Note that £2,2+3 < 0, and the result of plotting the entanglement factor Z = Y R /Y G — 
1 as a function of a(N) is presented in Figure 4. We see that the enhancement factor, 
which is fairly constant above a = 0.2008, begins to increase dramatically below this value 
of a = 0.2008. In general, there is an inverse relation between N and N e . The precise 
relation appears to be unknown, and may in fact depend on the details of production of 
the polymer. But we know that a must have a weak dependence on N e , so that the tensile 
modulus Yr has a predominantly inverse linear dependence on N e through the factor c e (as 
observed experimentallyS), with a multiplicative amplification factor having a much weaker 
dependence on N e . 

For values of a much less than 0.2008, the approximations utilized in our calculations do 
not hold (see section III), and we see a divergence in the enhancement factor. As a — > 0, 
extremely short length scales are encountered, and fluctuations become more important so 
that higher order diagrams will need to be considered. Our goal in this paper was mainly 
to tackle the transition to a state of entanglement, and this aim has been achieved. 



We postulated an extension of the Cahn-Hilliard functional to describe entanglement in 
polymers. We extended the Cahn-Hilliard functional with two terms. One is an attractive 
nonlocal term which describes the effect of entanglement, and the other a local repulsive 
term indicative of excluded volume interactions. We developed a simple model to connect 
the parameters of our theory with the parameters of the underlying chains. We showed in 
Appendix A how the extended functional can be derived using notions from gaue theory. 
Using field theoretic techniques to go beyond the Gaussian approximation, we showed that 
the onset of entanglement is a crossover phenomenon, signalled by the effective diffusion 
constant going to zero. A reasonable estimate for the critical partial chain concentration at 
which this crossover occurs was obtained. Finally, we also computed the tensile modulus 
of the system, which showed a dramatic increase below the critical value of N e discussed 
earlier. Limitations of the calculations were also discussed. 



I would like to acknowledge discussions with Brian Kendrick regarding the gauge theoretic 
formulation of the problem. The work described in this paper was done under the LDRD-CD 
program on polymer aging at the Los Alamos National Laboratory. 



In this appendix, we shall show how to derive Eqn.f|, along with the choice of parameters 
made in section II, using notions from gauge theory. Let us start with: 



V. CONCLUSIONS. 
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(A.l) 
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where the variables have been defined in section II. We will use this form as our starting 
point to generate a more complete energy functional using gauge invariance. 

From Eqn. [A.l| we see that Uq is invariant under global translations of c(s), i.e. under 



c(s) — > c(s) + h where h is a constant. And the appropriate group to consider is TV The 
physical origin of this group can be traced back to the fact that the quadratic (positive, semi- 
definite) form of the energy density is dictated by expanding the internal energy around a 
minimum, in a Landau-like fashion. Physicality of Ti transformations demands that c'+c e > 
0, where d denotes the deviation of the specie concentration from its average, and that 
number conservation is guaranteed. These physical constraints will be incorporated into the 
evaluation of the partition function, after gauging U , in a manner similar to applying gauge 
constraints in QFT. 

Our physical motivation for seeking local gauge invariance under Ti is the same as 
that of Yang and Millsil, and in quantum electrodynamics (QED), where one observes 
the invariance of the noninteracting Lagrangian under certain global transformations. One 
then demands covariance of the theory when these symmetry operations are local i.e., when 
the transformations are space-time dependent. A reason for this, as given by Yang and 
Mills, is that one can now freely interchange between the fields as one moves through space 
and time, while leaving the physics covariant. It is important to note that gauge theory 
in QFT is not a result of the fact that the phase of the field is not measurable. In fact, 
Aharonov and Bohm showed many years ago that the phase in quantum mechanics is indeed 
observable. Local transformations under Ti generate concentration fluctuations which arise 
from entanglement. 

Following Yang and Mills,@ local gauge invariance of u$ under Ti motivates us to define 
new fields b, which have invariance properties appropriate to TV We define a covariant 
derivative — > (^- + qrbi), where r = is the generator of Ti, q is a 'charge', or 
equivalently, a coupling constant, and the 6-fields are analogs of the magnetic vector potential 
in electrodynamics. In our previous theory of self-assembly0, gauge fields arose from the 
underlying covalent bonds between the two species in the system. In the present case, where 
we wish to describe entanglement, the gauge fields are to be thought of as arising purely from 
statistical considerations alone. On the other hand, chemical cross-linking would provide a 
physical origin for the 6-fields in vulcanized systems. The energy functional for the 6-fields 
is defined a la Yang and Mills, via the minimal prescription. With this, our original internal 
energy density is transformed into: 

(3u — > f3u = f3u + I3ui n t + (3u Y m (A. 2) 

where u int refers to the interaction energy density, and uym is the energy density as- 
sociated with the Yang-Mills 6-fields alone. Equivalently, we may define the total energy 
functionals associated with these energy densities: 

(3U ^(3U = f3U + (3U int + PU YM , 

where 

(3u int = J(c) • 6(s) + / 6(s) • 6(s) (A.3) 

with 
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J(c) = gqVc (A.4) 

(A.5) 




We need one more definition for completeness: 

^«=(i) [wr<£)wri£) s T (A - 6) 

This equation can be cast into the following form: 

(3u YM = - b t V\ (A.7) 

Eqn.|A.7| is obtained via an integration by parts, in the transverse gauge. Since we 



are dealing with an Abelian gauge theory, it is permissible to insert this transverse gauge 
manually, without resorting to the formal machinery of Faddeev and Popov. 
Again using the transverse gauge and integrating by parts, it is clear that: 

J d 3 sVc{s) ■ b{s) = - J d 3 s c{s) [V ■ b{s)} = = i J d 3 sVc{s) ■ b{s) (A.8) 

It is this crucial identity which allows us to get the precise form for Eqn.|], which we 
motivated in section II using an intuitive approach. It is the nature of the Ti group which 
permits this manipulation to go through successfully. In our earlier papertO, where we used 
the SO (2) group, such a manipulation would not have availed us any advantage. 

Note that we are utilizing a non-relativistic version of the Yang-Mills procedure, since we 
are only concerned with time-independent problems. Furthermore, since we are concerned 
with translations in Ti, there is only a single generator to contend with, so that the resulting 
functional is only quadratic and not quartic in the 6-fields. 

It is important to emphasize that the usual application of the Yang-Mills procedure in 
QFT implies the existence of fundamental interactions. In our case, we are applying the 
principle of local gauge invariance at the mesoscale. Consequently, we do not expect to 
discover any new fundamental interactions by using gauge invariance. Rather, we interpret 
the new 6-fields as yielding correlations between the concentration fields. These correlations 
could also be thought of as effective interactions, which arise at the mesoscale from the 
underlying electrostatic interactions between molecules. 

The partition function we need to evaluate is now: 

Q'= [Vc6 (c) J] Vb k exp - (3{U + U mt + U YM ) (A.9) 

fc=l,3 



Equation |A.9| is a functional integral, where the step functions denoted by 9 imply that 



we must restrict integration to positive semi-definite values of the fields. 

Since the 6-fields appear only quadratically in the above functional, it is straightforward 
to integrate over them, and obtain an effective internal energy functional involving only c, 
upon using Eqn. [A.8| . The result is: 
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(3U cS = (3U + (3AU eS = f3U + ^J d 3 s J dV J,(c(s)) ( -73 Tya ) J *( c ( s ')) ( A -!0) 



Note that in doing so, we have ignored an overall trivial normalization constant that 
appears in the evaluation of the partition function Q'. This is permissible, as this factor 
cancels during the evaluation of averages of observable quantities. 

To reveal the physics in this effective functional, we perform some straightforward algebra 
to write our result as: 



(a 2 \ r , fa 4 \ r „ r , , exp(- \2a 2 /g\s - s'|) , 
#7* = #/„ + ( y j / d 3 . c(s)c(s) - ( |-J y d 3 s J d 3 s' c(s) - V ^ c(s') 

(A.ll) 

where a 2 = g 2 q 2 /2 (we shall use units in which q=l). Note that U e s is quadratic, 
the generator of Ti making sure that higher order terms do not appear in our functional. 
Equation |A.11| is one of the main results of our paper, and provides a deeper motivation for 
the model developed in section II on intuitive grounds. As discussed in section II, the form 
of Equation |A.11| guarantees number conservation. Equation | A . 1 1| shows that entanglement 



may be understood in the context of a mesoscopic gauge theory. Note that the two terms we 
just discussed have signs opposite those of corresponding terms in theories of self-assembly0. 
We thus see that using Ti instead of SO (2) in the previous theory0 has led to a qualitatively 
different theory. 
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FIGURES 



FIG. 1. (a) is a pictorial representation of the cubic term in A. Each leg corresponds to a 
factor of c, the field. The intersection of the three legs symbolizes a factor of 7 = 1/6, the coupling 
constant, (b) is a pictorial representation of the quartic term in A. A factor of —1/12 is to be 
inserted at the intersection. 

FIG. 2. (a) represents the tadpole diagram which is crucial in our calculations, (b) represents 
the setting sun diagram. Both (a) and (b) are second order contributions to the correlation function 
coming from the cubic interaction term, the first order corrections being null. 

FIG. 3. This figure represents 1-loop (bubble) contribution from the quartic interaction term 
in A. 

FIG. 4. This is a plot of Z = Yr/Yq — 1 as a function of a. Notice that the factor is virtually 
constant above a = 0.22, followed by a dramatic increase below this value of a. Approximations 
employed in the calculation cause the entanglement factor to diverge at a ~ 0.01. Remember, 
decreasing a corresponds to increasing entanglement. 
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(b) 



Chitanvis, Fig. |3] 
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